Third minima in actinides - do they exist? 
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We study the existence of third, hyperdeformed minima in a number of even-even Th, U and Pu 
nuclei using the Woods-Saxon microscopic-macroscopic model that very well reproduces first and 
second minima and fission barriers in actinides. Deep (3^4 MeV) minima found previously by 
Cwiok et al. are found spurious after sufficiently general shapes are included. Shallow third wells 



may exist in 



Th, with Ilird barriers < 200 and 330 keV (respectively). Thus, a problem of 



qualitative discrepancy between microscopic-macroscopic and self consistent predictions is resolved. 
Now, an understanding of experimental results on the apparent third minima in uranium becomes 
an issue. 
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The topic of third minima in actinides, supposedly 
more deformed than the superdeformed (SD) ones, and 
for this reason sometimes called "hyperdeformed", is sur- 
rounded by some uncertainty. Their existence is inferred 
from an analysis of the observed transmission resonances 
in the prompt fission probability in (n,f), (d,pf), (7,f) and 
(e,f) reactions. Using one-dimensional models for tunnel- 
ing one tries to fit the data by considering various pos- 
siblities of rotational or vibrational states, either in the 
second or the third well, and choosing unknown barrier 
parameters. The best fit gives energy at the third barrier 
and minimum, and possibly the moment of inertia in the 
Ilird well. Additional check of the hypothesis may be 
obtained by measuring fission fragment angular distribu- 
tion. First interpretation of the third minima in thorium 
j? was supported by later studies [SHT; in 230,23i,232,233rpj^ 
which found them rather shallow ~300 keV. Recent ex- 
periments of the Debrecen-Munich group report third 
minima in a number of nuclei (232,234,236]q 232pg^^ [SLTTP. 
Based on resonance energies inbetween the first and the 
(higher) second barrier and their small widths, they infer 
a third barrier at least as high as the second one and deep 
third minima. There are not any direct measurements of 
the third minima so some reservations concerning the in- 
terpretation would be in order. 

Theoretical predictions of the third minima in ac- 
tinides present somewhat confusing picture. While the 
presence of the second minima at deformation close to 
j32 = 0.6 is a common feature of both macroscopic- 
microscopic (since Strutinsky [T) and self-consistent cal- 
culations, the third minima appear or disappear from 
calculation to calculation. The first calculated Ilird min- 
ima p^il3j were predicted shallow, 0.5 to 1.5 MeV, with 
the third barrier higher than the second one in thorium. 
Bengtsson et al. |T4] had very shallow minima (0.5 MeV 
or less), with the lower third barrier. Qualitatively dif- 
ferent results, with prominent 3-4 MeV deep, and often 
multiple, third minima, also in U and heavier nuclei, were 
presented by S. Cwiok et al. p!5l4T7] . These calculations 
used a richer set of axially symmetric deformations than 
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the previous ones, up to the multipolarity A = 7. Results 
of the Hartree-Fock calculations indicate that presence 
or absence of third minima depends very much on the 
effective interaction. However, minima in U and Pu, if 
they appear at all, are rather shallow and nearly degen- 
erate with, or even lower than, the ground states [2T] . 
Therefore, they do not fit Debrecen-Munich results. On 
the other hand, calculations with the Gogny force [22] 
and SkM* [23 (Hartree-Fock plus BCS) predict shallow 
minima in some Th isotopes, compatible with Th data 
[5 , while the other work [24 provides shallow minima at 
the top of the first barrier, not of the type suggested by 
the experimental papers. 

The aim of this work was to check the existence of 
the third minima in actinides within the Woods-Saxon 
microscopic-macroscopic model used in [151 US] by ex- 
tending the number of included deformations. An im- 
pulse to this came from first repeating old calculations 
and realizing that they produce HIrd minima also in 
many heavier actinides, in which they were not found 
experimentally. Moreover, the calculated minima with 
larger octupole deformations have quadrupole moments 
(calculated assuming the sharp-surface distribution and 
the radius parameter 1.16 fm) Q ~ 170 b, disturbingly 
close to the scission region. The contemplated possibility 
was that some Ilird minima in |T5l [16] are just interme- 
diate configurations on the scission path, whose energy 
was calculated erroneously because of limitations of the 
admitted class of shapes. Method Third minima were 
searched by exploring energy surfaces in a region of de- 
formations beyond the second minima. The energy is cal- 
culated within the microscopic-macroscopic method with 
a deformed Woods-Saxon potential. Nuclear shapes are 
defined in terms of the nuclear surface [27 

R{0,^) = c({/3})i?o(l + ^/3AlAo(^,¥') + 

A>1 

A>l,/i>0,even 

E /3x^Y^^s{0,^)), (1) 

A>l,/i>0,odd 

where c{{/3}) is the volume fixing factor, the real- valued 
spherical harmonics Yx^c^ with /i > even and Ya/^s, 



with /x > odd, are defined in terms of the usual ones as: 

Iapc = (iAM+rA-M)/\/2andyAM, = -i(yAM+yA-M)/\/2. 

Such shapes have at least one symmetry plane y — z. Al- 
though nonaxiality modifies energy landscapes in the re- 
gion of the third barrier, the main effect comes from a 
proper treatment of the axially symmetric shapes. There- 
fore, we explore here the first part of the Eq. (1). 

Two models were used for the macroscopic part: the 
Yukawa plus exponential [25^, used by us up to now in 
many studies, and recently given version of the Liquid 
Drop [2F with the mean-curvature term included in the 
deformation-dependent part of energy. The reason for 
using the latter was to check whether double third min- 
ima in ^15HT7] are not related to some singularity of the 
used Yukawa-plus-exponential macroscopic energy. 

The s.p. potential parameters used in the present 
work, as well as the way the shell- and pairing correc- 
tions are the same as in a number of our previous stud- 
ies, e.g. in p!8]42Q] . Let us emphasize that the used 
model very well reproduces first [18 , and second barriers 
[19 and second minima [20 in actinide nuclei. Removig 
shape limitation A general difficulty in studying very de- 
formed nuclear shapes lies in their effective parametriza- 
tion. The Eq. (IT]), natural for small and moderate defor- 
mations, becomes problematic for this purpose. Rather 
large (more than six) number of deformations is necessary 
for a sufficient accuracy. Their sizes may be misleading, 
for example, there are sets {Px} corresponding to huge 
elongations with quite moderate ^2- Such is the case of 
two kinds of Ilird minima in actinides shown in p!5UT7] . 
With similar j32 ^ 0.85 they might be thought close in 
elongation. But they are not: they correspond to differ- 
ent quadrupole moments Q ?^ 120 and 170 b. The energy 
mapping in many-dimensional space becomes a problem. 
A reduction of dimension via the minimization over some 
of the deformations often leads to an energy surface com- 
posed from disconnected patches, corresponding to mul- 
tiple minima in the auxiliary (those minimized over) di- 
mensions. If a smooth look of such a surface hides the 
jumps between minima, it suggests a false picture of the 
barrier. 

It may be for these difficulties that an important shape 
limitation of Eq. (IT]) did not attract much attention. The 
dipole deformation /3i is omitted there, as corresponding 
to a shift of the origin of coordinates which leaves energy 
(always calculated in the center of mass frame) invariant. 
However, this is true only for weakly deformed shapes. 
For large elongations, Pi acquires a meaning of a real 
shape variable. 

The parametrization ([T]) is not unique, as the same 
shape has many expansions corresponding to various 
shifts of the origin of coordinates. Hence, one could think 
that j3i may be always transformed out. However, the 
shift invariance is lost beacause a) the expansion (IT]) is 
always truncated, b) for well constricted shapes, large 
shifts lead outside the formula (IT]) making the radius mul- 
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FIG. 1: Energy along a sequence of smoothly elongating 
shapes parametrized by /3i — /^s, beginning at the apparent 
third minimum with larger octupole deformation /^s ~ 0.6 in 
^^'^U. Quadrupole moments calculated for the sharp-surface 
distribution and the radius parameter tq — 1.16 fm. De- 
picted shapes correspond to: Q = 171 (the apparent Ilird 
minimum), 189, 205 and 220 b. 



tivalued. In practice, admitting nonvanishing /3i in (IT]), 
R{e, ^) = c{{l3})Ro{l + Yl /^Ano(^)), (2) 

A=l 

turns out to be an effective way of including an im- 
portant class of smoothly constricted elongated shapes, 
which are unreachable when using a reasonably truncated 
expansion ([T]) without it. One should notice that the 
parametrization ([2| is redundant: the remaining invari- 
ance to small shifts of the origin of coordinates means 
that some shapes have many parametrizations with dif- 
ferent /3i. 

Disappearence of spurious minima We considered even- 
even nuclei 230,232Th, 232,234,236,238u and ^^Op^. Axially 
symmetric shapes were assumed. Two sets of energy cal- 
culations were done: A) maps (/32, Ps) with minimization 
over P^-Ps (7D) and B) in a hypercube Pi-Pe (6D) from 
which various 2D maps could be generated. The grid 
points were (with step in all directions = 0.05): —0.35-0 
in pi, 0.55-1.5 in p2, 0-0.35 in ps, -0.10-0.35 in p^, -0.2- 
0.2 in /^s, —0.15-0.15 in Pq. To improve limits on barrier 
heights we calculated energies along chosen continuous 
paths in 8D space Pi-Ps- 

Results of A) essentially reproduce landscapes in |T5l 
16 . Beyond the axially- and reflection-symmetric sec- 
ond minima, the reflection-asymmetry lowers the second 
saddles. Even more (shell correction of ~10 MeV) it 
lowers energy for still larger elongations. Double Ilird 
minima with Ps ~ 0.3 and 0.6 are obtained with the 
Yukawa-plus-exponential macroscopic energy; the less 
mass-asymmetric ones are deeper in Th and higher in 
U and Pu. With the LSD macroscopic energy, formula 
|26J minima with larger Ps (quadrupole moments) vanish 





and the Ilird minima become single. In both versions 
the third barriers are quite high, > 2.5 MeV above the 
Ilird minimum, with the highest ones in thorium. 

The landscape modification obtained by including 
dipole deformation /3i in ^ is decisive for the picture 
of the Ilird minima. The apparent minima with large 
Ps and quadrupole moments disappear completely in the 
studied nuclei. One can find continuous 8D paths start- 
ing at the supposed Ilird minimum and leading to scis- 
sion, along which energy decreases gradually. An exam- 
ple is shown in Fig. 1 for ^^^U. 

The existence of minima with smaller octupole defor- 
mation (and quadrupole moments) requires a more de- 
tailed study. Finding the barrier in many-dimensional 
space requires hypercube calculations. The maps (/32, Ps) 
generated by minimization over P^-Pq from the 6D hyper- 
cube calculations B) are shown in Fig. 2 for 232]j ^^ fixed 
/3i = and —0.2. The large effect of included new shapes 
may be seen as a lowering of the /3i = third barrier by 
Pi = —0.2. Its size of 2-4 MeV is common in all stud- 
ied nuclei, leaving some small barriers in 230,232rpj^ q^^ 
can note that the heights of the second saddles are prac- 
tically not changed by Pi. Hence, in the evaluation of 
secondary fission barriers in actinides (see e.g [19 ) this 
parameter is not necessarily needed. The results B) in- 
dicate a large effect of new shapes, but are not precise 
enough for a detailed barrier prediction. Our tests show 
that deformations Pi-Ps are needed for fixing minimima 
up to 100-150 keV, a greater precision requires at least 
Pq. As 8D or 9D hypercube calculations are very time- 
consuming, we find the upper limit on barriers selecting 
continuous paths in 8D Pi-Ps space, from the supposed 



Ilird minima towards scission. Energies vs quadrupole 
moments along such paths are shown in Fig. 3 for 232^-^]^ 
and ^^^U. As can be seen, the barrier vanishes in ura- 
nium and must be smaller than 330 keV in ^^^Th. The 
only other nonzero upper limit on the Ilird barrier of 200 
keV we find in ^^^Th. As by increasing the number of 
deformations one can only lower saddle points, the third 
barriers could rise only by lowering the supposed Ilird 
minima. But from our tests this effect is of the order of 
100 keV, so, within the used model, the barrier in ^32^^]^ 
should not be higher than ^ 450 keV. The results of this 
work may be summarized as follows: 

(i) With the shape parametrization ([l]), modifications 
of energy landscapes by including the dipole deforma- 
tion Pi become important beyond the Ilnd barrier in ac- 
tinides. 

(ii) They remove minima with larger octpole defor- 
mations found in [151 HI]- They also remove min- 
ima with smaller octupole deformations, perhaps, except 
230,232 rpj^^ where only a very shallow minima can exist. 
Time-consuming, 9D calculations are required for a pre- 
cise (up to less than 100 keV) determination of the Ilird 
barriers. 

(iii) The Ilird barriers predicted by the Woods-Saxon 
model, if they exist at all, do not exceed a few hundreds 
keV which is in agreement with selfconsistent mean-field 
results with realistic forces. They do not correspond to 
a type suggested originally in [151 US] which served a 
physical interpretation of some experimental data. 

It remains to understand the observed resonances in 
fission probability and their structure. On the theory 
side, with the mean-field barriers so shallow, one should 



check collective (beyond mean-field) corrections. One 
could also check barriers at higher spins. 




FIG. 2: Energy maps for ^^^U from the 6D /3i — /^e calculation, 
minimized over (3^^ Pb] upper panel - /3i =0, middle panel - 
/3i = —0.2. The lower panel gives the difference Eupper — 
Emiddiej SO the lowcriug of the third barrier by /3i = —0.2 can 
be read from it as a positive number. 
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FIG. 3: Energy along a sequence of continuously elongating 
shapes parametrized by /3i — /^s, beginning at the apparent 
third minima with smaller octupole deformation /Ss ~ 0.3 in 
^^^Th and ^^^U. Quadrupole moments calculated as in Fig.l. 
Depicted shapes correspond to states in ^^^Th: Q=117 (ap- 
parent Ilird minimum), 137 and 155 b. 



